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Abstract. We present intermediate- and low-resolution optical spectra around Ha and Li I A6708A for a sample 
of 25 low mass stars and 2 brown dwarfs with confirmed membership in the pre-main sequence stellar a Orionis 
cluster. Our observations are intended to investigate the age of the cluster. The spectral types derived for our target 
sample are found to be in the range K6-M8.5, which corresponds to a mass interval of roughly 1.2-0.02 Mq on the 
basis of state-of-the-art evolutionary models. Radial velocities (except for one object) are found to be consistent 
with membership in the Orion complex. All cluster members show considerable Ha emission and the Li I resonance 
doublet in absorption, which is typical of very young ages. We find that our pseudo-equivalent widths of Ha and 
Li I (measured relative to the observed local pseudo-continuum formed by molecular absorptions) appear rather 
dispersed (and intense in the case of Ha) for objects cooler than M3.5 spectral class, occurring at the approximate 
mass where low mass stars are expected to become fully convective. The least massive brown dwarf in our sample, 
S Ori45 (M8.5, ~0.02 Mq), displays variable Ha emission and a radial velocity that differs from the cluster mean 
velocity. Tentative detection of forbidden lines in emission indicates that this brown dwarf may be accreting mass 
from a surrounding disk. We also present recent computations of Li I A6708 A curves of growth for low gravities 
and for the temperature interval (about 4000-2600 K) of our sample. The comparison of our observations to these 
computations allows us to infer that no lithium depletion has yet taken place in a Orionis, and that the observed 
pseudo-equivalent widths are consistent with a cluster initial lithium abundance close to the cosmic value. Hence, 
the upper limit to the a Orionis cluster age can be set at 8 Myr, with a most likely value around 2-4 Myr. 

Key words, circumstellar matter — stars: abundances — stars: evolution — stars: late-type — stars: low mass, 
brown dwarfs — stars: pre-main sequence — open clusters and associations: a Orionis 
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1. Introduction 

Deep photometric and spectroscopic searches in various 
nearby star-forming regions and young open clusters have 
revealed populations of very low mass stars (<0.3M Q ), 
brown dwarfs (see Basri |200C for a review) with masses 
below the hydrogen-burning mass limit (^0.075 M Q ) and 



plane tary-m ass objects (Najita, Ti ede, fc Carr 2000; Lucas 
et al. 2001 ; Zapatero Osorio et al. |2000|) smaller than the 
deute rium-burning threshold at 0.013 Mq (Saumon et al. 
1996 ) . Age is one of the most relevant parameters for their 
study and characterization. Traditionally, mass estimates 
of cluster members rely on evolutionary models, which 
predict luminosities and effective temperatures (hereafter 
T e ff) as a function of time and mass. This procedure, 
namely isochrone fitting to observational photometric di- 
agrams, is highly age dependent for the smallest objects 
because t hey co old down very quic kly (e.g., D'Antona & 
Mazzitelli |l994| ; Burrows et al. |1997| ; Chabrier et al. |2000D . 
To better constrain masses, it is desirable to date clusters 
with higher accuracy. 

Lithium absorption features in optical spectra can be 
used as a tracer of the stellar internal structure. In ad- 
dition, lithium is valuable to assess the age of stars in 
clusters. Dating young clusters and field objects based on 
lithium analysis is a procedure nearly independent of dis- 
tance and reddening, in marked contrast to the isochrone- 
fitting technique. Moreover, lithium isochrones do not sig- 
nificantly depend on metallicity (DAntona 200C ), render- 
ing the lithium dating technique very powerful. Pre-main 
sequence stars with large convective regions burn lithium 
efficiently on short time scales (see Pinsonneault 1997 for a 
review) as soon as the temperature at the base of the con- 
vective zone becomes hot enough to undergo the nuclear 
reaction 7 Li + p — > 4 He + a. Stars smaller than the Sun re- 
quire only about 10-15 Myr to deplete this element by one 
order of magnitude, and all M-type stars are observed to 
have destroyed their lithium at ages around 20-40 Myr 
Pinsonneault, Kawaler 
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Fig. 1. Optical color-magnitude diagram of our a Orionis 
targets (large symbols) . Open circles indicate sources with 
Ha emissions larger than pEW=10A. Open triangles 
stand for cluster members with forbidden emission lines. 



S Ori objects (small crosses) from Bejar et al. (1999) 
are included in the figure for a better delineation of the 
a Orionis cluster photometric sequence. Overplotted onto 
the data are t he no- dust 5 Myr isochrone (dashed line) of 
Baraffe et al. ( 1998 ) and the 3 Myr isochrone (solid line) 
of D'Antona & Mazzitelli ( 1997 ), for which masses in solar 
units have been labelled. Spectral types as a function of 
the (R — I) color are also provided. 
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(e.g. 

& Mazzitelli |1994 |1997| Baraffe et al 
lithium detections in fully convective objects near and be- 
low the substellar limit have been successfully used to con- 
strain the a ges of clusters like the Plei ades (Basri, Marcy 
& Graham |l996| M artin et al. |l998[ Stauff er, Sc hultz, 
& Kirkpatrick |l998|) , aPersei (Stauf fer et al. |l999|) , and 
IC2391 (Barrado y Navascues et al. |l999| ). Lithium dat- 
ing, which is fundamentally a nuclear age calibrator, can 
be considered reliable even though some uncertainties (ro- 
tation, activity, mixing processes) may affect theoretical 
calculations. 

Recently, various groups have investigated the stel- 
lar and substellar populations around the bright, massive 
and multiple 09.5V-type star a Orionis, which gives its 
name to the association (Walter et al 



Demarque 1990: D'Antona 

1998j ). Furthermore, tween lMyr and 7 Myr (Blaauw |l96 



Walter, Wolk, & Sherry 1998; Bejar, Zapatero Osorio, k, 
1999| ; Zapatero Osorio et al. |1999| , p000| ; Bejar et 



1994: Wolk 1996 



California Institute of Technology, the University of California 
and the National Aeronautics and Space Administration (the 
Observatory was made possible by the generous financial sup- 
port of the W. M. Keck Foundation) . 



al. p001[) . These authors have adopted a cluster age be- 

Warren & Hesser 

1978| ; Brown, de Geus, & de Zeeuw |1994T) . This is the 
age interval estimated for the 09.5V-type star based on 
its physical properties, evolutionary stage (still burning 
hydrogen on the main sequenc e) and membership in the 
Orion OBlb subgroup (Blaauw 1991 ). Other properties of 
the a Orionis cluster, e.g., distance (352 pc) and redden- 
ing (Ay < 0.5mag), are discussed in Bejar et al. ( 2001 ). 
Here we examine low mass stars and brown dwarfs with 
confirmed membership to determine the most likely age 
of the cluster. We report on observations of intcrmediate- 
and low-resolution optical spectroscopy in Sects. ^ and El 
A discussion and main conclusions are given in Sects. Iq 
and |1, respectively. 



2. Sample selection 

Our list of 28 targets (12.3 < / <19.6, T cff ~ 4200-2400 K) 
comprises a Orionis solar-mass and low mass stars, and 
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Table 1. Log of the observations (ordered by observing date). 
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Table 2. Instrumental setups of different campaigns. 



Run 
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brown dwarfs selected from different surveys (see Table 0). 
All have been identified as genuine cluster members us- 
ing various techniques. Stars labelled with "4771" and 
"r" were first identified by Wolk ( fl996| ). VRIJHK pho- 
tometry, spectroscopy, and strong X-ray detections (in 
many cases) are available. Wolk (1996) provided equiv- 
alent widths of the Lil resonance doublet for a few of 
these stars. However, his spectroscopic data of relatively 
faint sources have poor signal-to-noise (S/N) ratios, which 
severely affects the measurements. We decided to re- 
observe them to improve the quality of the spectra. S Ori 
targets (IAU nomenclature) have been selected from the 



RIJ survey of Bejar (2001). They nicely fit in the cluster 
optical-infrared sequence. Albeit we lack previous spectro- 
scopic data for them, the spectra presented here confirm 
them as very active, young sources, and therefore, they 
have to be cluster members. The two brown dwarfs in our 
sample, S Ori 27 and 45, have been taken from Bejar et al. 
( 1999| ), where they are discussed at length. 



Names, coordinates and / magnitudes of our sample 
are provided in Table [l]. Their location in the optical 
color-magnitude diagram is illustrated in Fig, [l , where 
RI ph otome try has bee n take n from Wolk ( 1996 ), Bejar 
et al. (1999) and Bejar (2001). Overplotted are the solar 
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Fig. 2. CAHA spectra. Data have been normalized to 
the counts at around 7050 A and shifted by 1 in the upper 
panel and by 1.5 in the lower panel for clarity. 



metallicity, and "no-dust" 3Myr and 5Myr isochroncs of 
D'Antona & Mazzitelli ( |1997| ) and Baraffe et al. ( |199SD , 
respectively. Masses as predicted by the models are in- 
dicated in the figure. Co mpari sons with other tracks are 
provided in Bejar et al. ( 1999 ). Our targets have masses 
ranging from 1.2 Af Q down to roughly 0.02 A/©. 
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Fig. 3. ORM spectra. Data have been normalized to the 
counts at around 6535 A and shifted by 1 for clarity. 

3. Observations and data reduction 

We acquired intermediate- to low-resolution optical spec- 
tra using the following telescopes: the 3.5-m telescope 
at Calar Alto (CAHA) Observatory in Almeria (Spain), 
the 2.5-m Isaac Newton telescope (INT) at the Roque de 
los Muchachos Observatory (ORM) on the island of La 
Palma (Spain), the 10-m Keck II telescope at Mauna Kea 
Observatory on the island of Hawaii (U.S.A.), and the 
Otto Struve 2.1-m telescope at McDonald Observatory 
in west Texas (U.S.A.). Observing campaigns, spectro- 
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Fig. 4. Keck spectra of two brown dwarfs in the a Orionis 
cluster. Data have been normalized to the counts at 
around 7050 A and shifted by 1.5 for clarity. The spec- 
trum of S Ori 45 has been smoothed with a boxcar of 5 
pixels to increase the S/N ratio of the data. 



fa «* 

CD 



i — i 

a> 

CM 



1-053849-0238 



6200 6400 6600 6800 
Wavelength (Angstroms) 

Fig. 5. McDonald spectra. Data have been normalized to 
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graphs attached to the Cassegrain focus of each tele- 
scope, detectors, gratings and slit widths used for col- 
lecting data are summarized in Table |[ The red arm 
of the TWIN instrument (CAHA) and the 235 mm 
camera of the IDS spectrograph (ORM) were chosen. 



Table [l] shows the journal of the observations, which in- 
cludes the nominal dispersions of the instrumental se- 
tups. No binning of the pixels along the spectral di- 
rection and the projection of the slits onto the de- 
tectors yielded spectral resolutions of 1.54 A (i?~4400, 
CAHA, first run), 4.82 A (i?~1600, CAHA, second 
run), 1.68 A (-R-3800, ORM), 2.89A (i?-2500, Keck), 
1.4 A (i?^4600, McDonald), and spatial resolutions as 
listed in Table ^. A binning of 8 pixels along the spatial 
direction was applied to the CCD at McDonald. Filters 
blocking the light blueward of 5000A were used at the 
CAHA and ORM telescopes. No order-blocking filter was 
used at the Keck II telescope; nevertheless, the two tar- 
gets observed are very red and the contribution of their 
blue light to the far-red optical spectrum is negligible 
(Martin et al. 1999| ). Weather conditions during the four 
runs (CAHA, ORM, Keck and McDonald) were spectro- 
scopic. The seeing in the visible was stable at around 1" at 
Keck, and 1.5"-2" at CAHA and ORM. The spatial reso- 
lution at McDonald was 2.72"/pix due to some technical 
problems related to the instrumentation. 

Raw images were reduced with standard procedures 
including bias subtraction and flat-fielding within NOAO 
irafQ We extracted object and sky spectra using the op- 
timal extraction algorithm available in the APEXTRACT 
package. A full wavelength solution from calibration lamps 
taken immediately after each target was applied to the 
spectra. The mis of the fourth-order polynomial fit to 
the wavelength calibration is typically 5-10% the nominal 
dispersion. To complete the data reduction, we corrected 
the extracted spectra for instrumental response using data 
of spectrophotometric standard stars (HD 19445, Feige34, 
G 191 B2B, BD+26 2606) obtained on the same nights and 
with the same instrumental configurations. These stars 
ha ve flux es available in the iraf environment (Massey et 
al. (19881) . 

The resulting spectra are depicted in Figs. ^]-|| They 
are ordered by increasingly late spectral type and shifted 
by a constant for clarity. The region around the Li I 
A6708A line is amplified in Figs. ||-[|. In Fig. [?| we have 
included the spectrum of the field M6-type spectroscopic 
standard star G1406 for a better comparison. 

4. Analysis and results 

4.1. Spectral types 

We inferred spectral types by comparing our target spec- 
tra to data of spectroscopic standard stars (G1820A, K5V; 
G1820B, K7V; G1338A, MOV; G1182, M0.5V; G1767A, 
M1V; G1767B, M2.5V; G1569A, M3V; G1873, M3.5; 
G1402, M4V; G1905, M5V; and G1406, M6V). The refer- 
ence spectra were obtained with similar instrumentations 
in previous campaigns. In addition, we observed several K- 



1 IRAF is distributed by National Optical Astronomy 
Observatory, which is operated by the Association of 



Universities for Research in Astronomy, Inc. 
with the National Science Foundation. 



under contract 
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and M-type standards with the McDonald telescope. For 
the spectral classification, we also used molecular indices 
that are based on the relative strengths of CaH and TiO 
bands (Kirkpatri ck, H enry & McCarthy 1991 ; Prosser, 
Stauffer & Kraft 1991), and the pseudocontinuu m PC 3 
index given in Martin, Rebolo & Zapatero Osorio ( 1996|) , 
which is valid for types later than M3. Our measurements, 
with an uncertainty of half a subclass, are provided in 
Table ||. The spectral t ypes o f S Ori 27 and 45 have been 
taken from Bejar et al. ( 1999 ). The final adopted spectral 
classes are in the range K6-M8.5. 

We note that our spectral classification relies on field 
dwarf objects with high gravities. The gravity of a Orionis 
cluster members is expected to be around log g — 4.0 (CGS 
units) accor ding to the evolutionary mode ls of Baraffe 
et al. ( |1998|) and DAntona & Mazzitelli ( |l994|) . Older 
K-type stellar counterparts in the field (^5Gyr) display 
similar gravities, but early-M and late-M stars have val- 
ues 0.5 dex and 1.0 dex larger, respectively. Cool giants are 
characterized by very low gravities (log g= 1.5-2, Bonnell 
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& Tell |1993[ van Belle |1999| ). Therefore, it is reasonable to 
base the spectral classification of young late-type objects 
on a scheme intermediate between that of dwarfs and that 



of giants. Luhman ( 1999 ) successfully applied this exercise 
to members of the young cluster IC 348, inferring that the 
spectral classification of objects like those of a Orionis can 
be obtained from dwarfs with an accuracy up to half a 
subclass. We have confirmed this by comparing the opti- 
cal spectrum of our M8.5 brown dwarf with brown dwarfs 
of identi cal ty pes in p O ph an d IC 348 (Luhman, Liebert 
& Rieke 1997 ; Luhman 1999 ). The three spectra overlap 
very nicely. We are confident that the spectral types given 
in Table [| are reliable within the quoted uncertainty. 

Since the spectral classification reflects effective tem- 
peratures, cluster members should lie along a defined 
sequence in magnitude vs spectral type diagrams. The 
a Orionis spectroscopic sequence is depicted in Fig. [l0[ 
where we have combined data presented here with data 
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a Note the drop of "S Ori" for some objects. 

6 Uncertainty of half a subclass. 

c Modified Julian date at the beginning of the exposure. 

d In emission. Whenever more than one spectrum available, the pEW has been measured over the combined data. 

e In absorption. Whenever more than one spectrum available, the pEW has been measured over the combined data. 



taken from Bejar et al. ( 1999| ), Barrado y Navascues et 
al. ( |2001a| ) and Martin et al. (|2001[ ). We note that the 
figure covers a wide range of masses: stars, brown dwarfs 
and planetary-mass objects. Free-floating low mass stars 
and isolated planetary-mass objects in the a Orionis clus- 
ter have luminosities in the /-band that differ by about 3 
orders of magnitude. Because substellar objects contract 
and fade very rapidly, such a difference becomes incredibly 
large at older ag es, e. g., 8 orders of magnitude at 100 Myr 
(Chabrier et al. Em]). 



4.2. Rotational velocities 

Given the poor velocity resolution of our spectra 
(68kms -1 , CAHA, first run; 176kms -1 , CAHA, second 
run; 78kms -1 , ORM; 120 kms -1 , Keck; and 65kms -1 , 
McDonald), we are able to detect extremely fast rota- 
tors. These are defined as objects with projected rota- 



tional velocities, vsini, larger than 55 kms - . The only 
case in Table || is 4771-1097, which was observed with 
the largest dispersion. This star (K6) is the most massive 
object in our sample (~0.9-l.2M©). We measured the ro- 
tational velocity by comparing its spectrum to a slowly ro- 
tating template selected from our sample. The spectrum 
of 4771-1075 (observed with the same instrumentation) 
has very sharp lines and a similar spectral type, as can 
be seen from Fig. ||. The full width at half-maximum of 
the atomic lines of this "reference" star indicates that its 
spectral broadening is mainly instrumental. Our proce- 
dure was to produce a set of artificial spectra spinned up 
to velocities of 75, 80, 95 and llOkms -1 . We then com- 
pared the observed spectrum of 4771-1097 to the synthetic 
rotational spectra and performed an analysis using the 
minimum squares technique. Spectral regions free of emis- 
sion lines and telluric absorptions were considered. We de- 
rived usini = 80±15 kms -1 . Wolk (1996) found clear evi- 
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dence of optical photometric variability in 4771-1097 and 
inferred a likely rotational period of ~lday. Combining 
these results with predicted radii for masses in the range 
0.9-1.2 Mq and ages between 3 and 5Myr (D'Antona & 
Mazzitelli |1994| ; Baraffe et al. |1998[ ), we conclude that 
4771-1097 is rotating with an inclination of i = 50°-90°. 

4.3. Radial velocities 

We computed radial velocities via Fourier cross- 
correlation of the target spectra with templates of sim- 
ilar spectral type. Our CAHA and ORM measurements 
were calibrated with the radial velocity standard star 
G114 (?;,. = 3.3 ±0.3 km s' 1 , Marcy & Benitz |1989| ; Marcy 
& Chen p^92| ), which was observed with the largest dis- 
persion at CAHA on 1998 Nov. 20. We used this star to 
derive the radial velocity of 4771-1051, and then corre- 
lated the rest of our CAHA targets of the same resolution 
against it. We also adopted 4771-1051 as the template for 
the ORM data. Whether this star has a variable radial ve- 



locity is unknown to us. Thus, our ORM radial velocities 
might be shifted by a certain amount. However, this is un- 
likely (at least within the error bars of the measurements) 
since there is another a Orionis member, r053833-0236, 
observed at CAHA and ORM. After correlating the ORM 
spectra of these two stars, we obtained for r053833-0236 
a heliocentric radial velocity similar to the one derived 
from the CAHA data. We note that the relative radial 
velocity of each target with respect 4771-1051 is reliable. 
We adopted r053833-0236 as the reference star for the 
low-resolution CAHA spectra, and the M8 field star vB 10 
was used as the template (v r = 35.3il.5kms" 1 , Tinney & 
Rcid 1998) for the Keck spectra. The spectrum of vB 10 
was taken from Martin et al. ( 1996| ). Our McDonald spec- 
tra were cross-correlated against the radial velocity stan- 
dard stars G1182 (v r = 32.4il.5kms" 1 , Jeffries |l995|) and 
Gl 87 3 (u r = 0.47±0.24kms" 1 , Marcy, Lindsay & Wilson 
1987), which were observed with the same instrumenta- 
tion and on the same nights. 
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Radial velocities, their uncertainties and the tem- 
plates used are provided in Table ||. We took special care 
in cross-correlating spectral windows (e.g. 6100-6800 A, 
8400-8800 A) that are not affected by telluric absorptions 
and that contain many photospheric lines. In addition, 
we considered only parts of the spectra free of emission 
lines. The error bars in the table point to a possible 1/4 
pixel uncertainty in the Fourier cross-correlation tech- 
nique (Martin et al. |l999t Lane et al. |200l[). We have 



Lane et al. 

checked this by cross-correlating the McDonald spectra 
against two reference stars. The spectrum of SOri45 is 
rather noisy, and the quoted error bar comes from the dis- 
persion observed after cross-correlating different spectral 
regions. The majority of our radial velocities are obtained 
to an accuracy of the order of lOkrns -1 . After discard- 
ing the largest and smallest radial velocity values from 
Table|(i.e., r053820-0237 and S Ori45, respectively), the 
mean heliocentric radial velocity of our a Orionis sam- 
ple is < v r > = 37.3 km s -1 with a standard deviation 
of 5.8 km s -1 . This is comparable to the systemic radial 
velocity of the cluster's central star, which has been de- 
termined to be in the range 27- 38 km s -1 (Bohannan & 
Garm any |1978| ; Garmany et al. |l980t Morrell & Levato 
199 1| ) . Additionally, these velocities (except for one, see 



Sect. H) are consistent with our sa mple b elonging to the 
Orion OB association (Alcala et al. [2000 ), and their dis- 
tribution is significantly different from that of field stars. 



4.4. Ha emission 

We derived Ha pseudo-equivalent widths via direct inte- 
gration of the line profile with the task SPLOT in iraf. We 
note that given the cool nature of our sample, equivalent 
widths in the optical are generally measured relative to 
the observed local pseudo-continuum formed by (mainly 
TiO) molecular absorptions (Pavlenko 1997| ). We will refer 
to these equivalent widths as "pseudo-equivalent widths" 
(pEWs). 

Because of the resolution of our observations, broad 
Ha lines appear blended with other nearby spectral fea- 
tures. The results of our measurements, given in Table |[ 
have been extracted by adopting the base of the line as the 
continuum. The error bars were obtained after integrating 
over the reasonable range of possible continua. Although 
this procedure does not give an absolute equivalent width, 
i.e., measured with respect the real continuum, it is com- 
monly used by various authors, and allows us to com- 
pare our values with those published in the literature. We 
note that all of our program objects show Ha in emission 
and that no significant Ha variability is found in any of 
them, except for r053833-0236 and S Ori45. We also note 
that the Ha emission of the fast rotator 4771-1097 is not 
stronger than that of other similar- type cluster members. 

Figure [n] shows the distribution of Ha pEWs as a 
function of spectral type. Effective temperatures are given 
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on the basis of the temperature — spectral-type relation- 
ships by Legg ett et al. ( 1996 ), Jones et al. ( 1995p and 
Bessell (1991). Masses as inferred from the 5Myr evo- 
lutionary isochrone of Baraffe et al. (1998) are also in- 
dicated in the figure. In general, there is a trend of in- 
creasing Ha emission for cooler spectral classes, i.e., for 
lower masses. This behavior has been observed in various 
young clusters, like the Pleiades and Hyades (Stauffer et 
al. |1994[ ), IC4665 (Prosser |l993| ) , aPersei (Prosser |l99§, 
and P raesepe (Barrado y Navascues, Stauffer & Randich 
199§| ). The relative increase of Ha in M-dwarfs may be 



(at least partially) explained by the drop of the flux con- 
tinuum and the larger TiO molecular absorptions in the 
optical as a consequence of cooler T e ffS. We note that, on 
average, Ha for a given spectral type is slightly larger in 
a Orionis than in other open clusters. This is very likely a 
direct consequence of the marked youth of a Orionis. 

In Fig. |ll| Ha emission appears very strong 
(pEWs > 20 A) and dispersed for late spectral classes 
(>M3.5), corresponding to masses below 0.25 M Q at the 
age of a Orionis. Various authors have found an appar- 
ent "turnover" in the distribution of Ha emission in the 
Pleiades (Stauffer et al. [1994 ; Hodgkin, Jameson & Steele 



1995) and aPersei (Zapatero Osorio et al. 1996). Pleiades 
and a Per stars with spectral types later than M3.5-M4 
show a lower level of emission than stars with warmer 
classes. The authors suggest that this turnover is due to 
the transition from radiative to convective cores. By in- 



specting D'Antona & Mazzitelli (|1994[ ) pre-main sequence 
evolutionary models, we find that this transition takes 
place at masses 0.3-0.2 M Q regardless of age. In a Orionis 
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we do not see a drop in the Ha emission of fully convec- 
tive objects, but an enhacement. The source of such large 
emission clearly diminishes by the age of the a Persei clus- 
ter (90Myr, Stauffer et al. 1999). However, the emission 
level of more massive stars remains with similar strengths. 

Three stars in our sample, namely 4771-41 (K7), 
SOriJ054001. 8-022133 (M4) and 4771-899 (K7), show 
profiles of Ha emission similar to those of classical T Tauri 
(CTT) stars, i.e., double peak structure and very broad 
lines spanning over ±300kms _1 from the line center. We 
illustrate in Fig. |l2| the region around Ha for these ob- 
jects. While the emission intensity is rather large in 4771- 
41 and SOriJ054001.8-022133 (pEWs above 45 A), it is 
moderate in 4771-899. 

We have calculated the Ha luminosity (Lhc<) for our 
sample as in Herbst & Miller (1989) and Hodgkin et 
al. ( |1995[ ). The ratio of Lh q to bolometric luminosity 
(Lna/iboi) is independent of the surface area and repre- 
sents the fraction of the total energy output in Ha. To de- 
rive Lboi we have used bolometric corrections provided by 
Monet et al. fll992j ) and Kenyon & Hartmann ( |1995[ ). The 
logarithmic values of Lua/L^i are listed in Table^; uncer- 
tainties take into account errors in photometry and in Ha 
pEWs. Figure [l3| shows the distribution of log (Lua/Lboi) 
with spectral type. For comparison purposes, we ha ve als o 
included the Pleiades mean values (Hodgkin et al. 1995 ). 
In the Pleiades, the Lua/Lboi ratio clearly increases to 
a maximum at around the M3 spectral type and then 
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[Oi] [Nil] He I [Sn] 

Object MJD" A6300A A6364A A6548A A6583A Ha A6678A A6716A A673lA 

(-51000) (A) (A) (A) (A) (A) (A) (A) (A) 



r053833-0236 


138.1237 


2.29±0.10 





75±0.05 


0.80±0.05 


2.50±0.10 


14.0±2.0 




<0.1 


1.15±0.08 


1.65±0.08 




138.1433 


2.00±0.10 





57±0.05 


0.88±0.05 


2.54±0.10 


14.1±2.0 





15±0.05 


0.85±0.05 


1.49±0.05 


J053949.3-022346 


138.1696 


1.80±0.50 




<0.2 


<1.0 


0.35±0.08 


42.0±6.0 


1 


18±0.08 


<0.15 


<0.15 


SOri45 6 


168.4644 






<1.5 


<1.5 


5.0±3.0 


33.0±9.0 




<1.5 


5.0±3.0 


3.56±2.0 


J054001. 8-022133 


204.0382 


1.27±0.10 





33±0.05 


0.26±0.05 


0.33±0.05 


47.0±9.0 


1 


05±0.08 


<0.15 


<0.15 




204.0531 


0.45±0.10 




<0.1 


0.10±0.05 


0.15±0.05 


46.0±9.0 





38±0.08 


<0.15 


<0.15 


4771-41 


205.9222 


0.96±0.10 




<0.1 


blended 


0.25±0.05 


53.2±9.0 





21±0.05 


0.16±0.08 


0.42±0.05 




205.9366 


1.00±0.10 





18±0.08 


<0.1 


0.30±0.05 


54.0±9.0 





42±0.07 


<0.1 


0.45±0.05 


r053840-0230 


206.0299 


<0.15 




<0.15 


0.42±0.05 


1.10±0.10 


6.5±0.6 




<0.1 


<0.15 


0.31±0.05 




206.0582 


0.55±0.05 





15±0.05 


0.54±0.05 


l.lOdhO.10 


6.9±0.6 




<0.1 


<0.15 


0.31±0.05 


r053849-0238 


515.2615 


0.36±0.10 




<0.1 


<0.1 


0.32±0.05 


2.6±0.3 




<0.1 


<0.1 


<0.1 


J053827.4-023504 


517.4297 


1.25±0.50 




<0.3 


<0.3 


<0.3 


21.2±3.0 





26±0.08 


<0.3 


<0.3 



a Modified Julian date at the beginning of the exposure. 

b Measures over the combined spectrum. Individual Ha pEWs were 20.0, 49.4 and 22.0±7.0 A, respectively. 



turns over. This is not observed in the a Orionis clus- 
ter, where cooler objects present larger Ha output fluxes 
than the older Pleiades spectral counterparts. Discarding 
a Orionis members with log (Ln a / Lhoi) > -3.2 dex, cluster 
data appear to display a flat distribution from late K to 
late M (i.e., no dependence on color and mass) at around 
log (Lua/Lboi) = -3.61 dex, with a standard deviation of 
0.18 dex. 



4.5. Other emission lines 

Our targets are pre-main sequence objects characterized 
by significant Ha emission and the presence of lithium 
in their atmospheres (see Sect. 4.6.2). All show proper- 



ties that resemble T Tauri stars. The nominal definition 
of weak-lined T Tauri (WTT) stars is usually based on 
Ha emission: pEWs smaller than 10 A for K and early-M 
stars (Herbig & Bell | 1988| ) and smaller than 20 A for later 
M- types (Martin 1998 ). This is accomplished by many of 
our objects. 

Some of our program targets display, however, other 
permitted (Hei A6678A) and forbidden ([Oi] A6300A, 
[Nil] A6548, A6583A, [Sn] A6716, A6731 A) emission lines. 
We have measured their pEWs; values are given in Table | 
as a function of Julian date. We note that some contam- 
ination from terrestrial night-sky emission lines may be 
expected in the measurements of the faintest sources. The 
objects of Table [| are plotted with different symbols in 
various figures of this paper, except for r053833-0236 (for 
this star we have used the "quiet" ORM data). The major- 
ity of the targets from Wolk (1996) are, in addition, clas- 
sified as strong X-ray emitters by this author. In contrast 
to the younger CTT stars, WTT objects are not accreting 
mass from disks. However, the presence of Hei and [Oi], 
[Nil], [Sll] emission lines is related to jets and outflows, 
which are typica l of C TT stars and accretion processes 
(Edwards et al. 1987; Hartigan, Edwards & Ghandour 



1995 ). These lines are generally detected in objects with 
strong Ha emissions (pEWs > 10 A, see Fig. |ll|). The co- 
existence of a Orionis members with properties of WTT 
and CTT stars is indeed indicative of ages of a few Myr. 
It may also indicate that small objects are accreting for 
longer periods than are more massive stars (Hillenbrand et 
al. 1998; Haisch, Lada & Lada 2001), provided that their 



strong Ha emissions are due to disk accretion. 

The star r053833-0236 shows strong Ha emission and 



[Nil] and [Sn] in two 
upper panel). However, 



noticeable forbidden lines of [O I 
consecutive CAHA spectra (Fig. |, 
its Ha intensity clearly decreased, and no other emis- 
sion lines were present in data collected on the following 
night (Fig. |, lower panel) or with the INT (Fig. |). The 
sources of this episodic flarelike event are not continuous 
in r053833-0236, probably indicating inhomogeneus mass 
infall onto the star surface. 

The case of the brown dwarf SOri45 (~O.O2M ) is 
particularly interesting and noteworthy. Albeit the detec- 
tion of [Nil] and [Sn] emission lines is affected by large 
uncertainties because of the modest quality of the Keck 
spectrum, this finding is very encouraging. It suggests that 
substellar objects, even those with very low masses, can 
sustain surrounding disks from which matter is accreted. 
Muzerolle et al. ( [2000 ) has recently reported on the evi- 
dence for disk accretion in a T Tauri object at the substel- 
lar limit. The presence of disks around brown dwarfs in the 
Trapezium cluster (~lMyr) has been proved by Muench 
et al. (2001). The emission lines observed in SOri45 indi- 
cate that "substellar" disks can last up to ages like those 
of the a Orionis cluster. It is also feasible that the probable 
binary nature of SOri45 (see Sect. ||) triggers the forma- 
tion of these emission lines. Nevertheless, further spectro- 
scopic data will be very valuable to confirm the presence 
of forbidden emission lines in S Ori45. The rapid Ha vari- 
ability of this brown dwarf is also remarkable. 
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4.6. Lil absorption 
4.6.1. Synthetic spectra 

We have computed theoretical optical spectra in the wave- 
length range 6680-6735 A around the Lil A6708A reso- 
nance doublet for gravity log g — 4.0 (CGS units) and for 
T off = 4000-26 00 K by running the WITA6 code described 
in Pavlenko ( 2000 ). This code is designed to opperate 
in the framework of classical approximations: local ther- 
modynamic equilibrium (LTE) , a plane-parallel geometry, 
neither sources nor drops of energy. The synthetical spec- 
tra have been obtained using the atmospheric structure of 
the NextGen models published in Hauschildt, Allard, & 
Baron ( 1999 ). We have adopted a microturbulent velocity 
value of v t = 2km s~ 1 , solar elemental abundances (Anders 



& Grevesse 1989 ), except for lithium, and solar isotopic 
ratios for titanium and oxygen atoms. The ionization- 
dissociation equilibria were solved for about 100 differ- 
ent species, where constants of chemical equilibrium were 
taken from Tsuji flL973p and Gurvitch et al. fll979| ). For 
the particular case of the TiO molecule, we have adopted 
a dissociation potential of £>o = 7.9eV and the molecu- 
lar line list of Plez (1998). The atomic line parameters 
have been taken from the VALD database (Piskunov et 
al. 1995), and the procedure for computing damping con- 
stants is discussed in Pavlenko et al. (1995) and Pavlenko 



(2001) 



Synthetic spectra were originally obtained with a step 
of 0.03 A in wavelength, and were later convolved with ap- 
propiate Gaussians to match a resolution of 1.68 A, which 
corresponds to the majority of our data. We have pro- 
duced a grid of theoretical spectra for nine different abun- 
dances of lithium [log A(Li) = 1.0, 1.3, 3.1, 3.4, re- 
ferred to the usual scale of log N (H) = 12] and seven values 
of T cff (4000, 3600, 2400, 3200, 3000, 2800 and 2600 K), 
covering the spectral sequence of our program tar- 
gets. Determinations of the meteoritic lithium abundance 



; Grevesse & Sauval |1998|) 
and 3.4. Extensive lithium 



(Nichiporuk & Moore |1974 
lie between logiV(Li) = 3.1 
studies performed in solar metallicity, intermediate-age 
clusters like the Pleiades (Soderblom e t al. 1993), a Per 
(Balachandran, L amber t & Stauffer 1996| ), 
(Jeffries & James |l999| ), NGC 2516 (Jeffries 
Thurston |l998|), and IC2602 and IC2391 (Randich et 



Blanco 1 
James & 



al. EOOlD , as well as in the Taurus star-forming region 
(Martin et al. 1994 ) show that non-depleted stars pre- 
serve an amount of lithium compatible with a logarith- 
mic abundance between 2.9 dex and 3.2 dex. We will adopt 
the mean value of log Ao(Li) = 3.1 as the cosmic "initial" 
lithium abundance. 

Figure [l4| depicts some of our theoretical spectra for 
different values of lithium abundance and surface temper- 
ature. The observed spectrum of S Ori 27 is compared to a 
few computations in Fig. [ll]. Optical spectra at these cool 
temperatures are clearly dominated by molecular absorp- 
tions of TiO. Only the core of the lithium line is observ- 
able, since the doublet wings are completely engulfed by 
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■ / v w: / \ 
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6700 6710 
Wavelength (Angstroms) 




6700 6710 6720 

Wavelength (Angstroms) 

Fig. 14. The upper panel shows theoretical spectra com- 
puted for Toff — 3400 K, log g — 4.0 and three different 
lithium abundances. The lower panel illustrates spectra 
for logg = 4.0, log A^(Li) = 3.1 and various temperatures. 
The spectral resolution is ~1.7A. The absorption feature 
centered at 6707.8 A is due to the atomic Lil resonance 
doublet, while the rest of the spectral features are mainly 
molecular TiO absorptions. 

Table 5. LTE Lii A6708A resonance doublet curves of 
growth: predicted pseudo-equivalent widths (pEWs). 











logiV(Li) 




(K) 


1.0 


1.6 


1.9 


2.5 


3.1 


3.4 


2600 


.357 


.444 


.479 


.552 


.617/.644* 


.656/.694* 


2800 


.346 


.440 


.475 


.551 


.623/.675* 


.669/.728* 


3000 


.312 


.404 


.441 


.522 


.596/.666* 


.652/. 741* 


3200 


.296 


.386 


.423 


.504 


.578/.639* 


.637/.729* 


3400 


.266 


.350 


.385 


.456 


.544/.634* 


.604/.720* 


3600 


.262 


.337 


.373 


.442 


.536/.566* 


.609/.665* 


4000 


.189 


.281 


.319 


.403 


.507/. 537* 


.588/.620* 



NOTES 
used log 
labelled 



— pEWs are given in A. In all computations we have 
g = 4.0 and solar metallicity, except for the columns 
with an asterisk, where we have used log g = 4.5. 



TiO lines (Pavlenko 1997). We have obtained the theoret- 
ical Lii A6708A pEWs via direct integration of the line 
profile over the spectral interval 6703.0-6710.8 A. Many of 
the lithium LTE curves of growth employed in this work 
are presented in Ta ble pi. Various auth ors (e. g., Magazzu, 
Rc bolo fc Pavlenk o |l992| ; Martin et al. |l994[ Pavlenko et 
al. 1995 ; Pavlenko 1998| ) have shown that the differences 
between LTE and non-LTE calculations for cool tempera- 
tures are negligible compared to uncertainties of pEW, T c g 
and gravity. Similarly, the effects of chromospheric activ- 
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Fig. 15. Synthetic spectra (T eff = 3000 K, log AT (Li) = 3.1) 
compared to the observed spectrum of the brown dwarf 
SOri27 (thick dotted line). Computed spectra have been 
degraded to the same resolution as the observations. The 
location of the Li I resonance doublet is indicated with a 
vertical dotted line. 



ity on the line formation are found to be of secondary im- 
portance (Pavlenko et al. 1995; Houdebine & Doyle 1995; 
Pavlenko 1998| ) and have not been included in our calcu- 
lations. The Li I resonance doublet appears to have very 
light dependence on the temperature structure of the outer 



layers (see also Stuik, Bruls & Rutten 1997). We find a 
rather poor agreement between the predicted Lil pEWs 
of Ta ble || and those provided in Pavlenko & Magazzii 
(|l996p . These authors' values are considerably larger be- 
cause they measured theoretical equivalent widths (note 
the drop of "pseudo" ) relative to the computed "real" con- 
tinuum, while we have determined pEWs relative to the 
computed pseudo-continuum formed by molecular absorp- 
tions. 



4.6.2. Observed spectra 

We have also obtained the Li I A6708 A pEWs from our ob- 
served spectra. To compensate for the different resolution 
of the data, the integration of the line profile has always 
been performed over the spectral range 6703.0-6710.8 A. 
Our measurements and their uncertainties are listed in 
Table |. Li I is detected in absorption in all of our pro- 
gram objects, except for S Ori J053914.5-022834 (M3.5). 
It might be a cluster non-member, but its optical spectrum 
is the noisiest amongst the McDonald data, and even the 
Cai line at 671 7 A lies barely undetected (see Fig. ^j. We 
impose a la upper limit of pEW = 0.44A by considering 
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Spectral type 

Fig. 16. Pseudo-equivalent widths of Lil A6708A as a 
function of spectral type. Symbols are as in Fig. 0. Note 
that the coolest cluster member of our sample (the brown 
dwarf S Ori 45) is not included in the figure for clar- 
ity. Overplotted onto the data are three LTE theoreti- 
cal curves of growth provided in this paper. Typical un- 
certainty in spectral type is half a subclass. The stellar- 
substellar borderline takes place at M5-M6 spectral type 
at the age of the cluster. Effective temperatures in Kelvin 
and masses in solar units are also given. 



the strongest possible feature in the region around the line 
and taking into account the S/N ratio and resolution of 
the spectrum. 

Lil pEWs are plotted against spectral type in Fig. [l6| 
SOriJ053914.5-022834 is excluded from the diagram. 
Overplotted onto the data are the theoretical pEWs 
for log g = 4.0 and two different lithium abundances: 
logiVo(Li) = 3.1 ("initial") and logTV(Li) = 1.9 (about one 
order of magnitude of destruction) . We have also included 
in the figure the "initial" curve of growth for a slightly 
larger gravity, logg = 4.5. The trend of the observations 
is nicely reproduced by the log-/Vo(Li) curves, implying 
that lithium is still preserved at the age of the a Orionis 
cluster. We will discuss this issue further in section |5.2| . 
We note the differences due to gravity in the Lil curves 
of growth. Although these differences are rather small 
(pEW<0.03A) for T cff >3700K, they become twice as 
large for cooler temperatures. Given the error bars of the 
observed Li I pEWs, we cannot easily discriminate between 
gravities. 

The scatter of the Lil pEWs is considerable for 
spectral types cooler than M3.5 (Fig. |l6|). The prob- 
lem of the lithium star-to-star dispersion occurring at 
T e ff < 5300 K has been widely discussed in the literature 
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Soderblom et al. 



199E; Randich et al. 



199$ Pallavicini et al. [199$ Russell 



1998; Barrado y Navascues et al. 



2001b). Nevertheless, this phenomenon still remains ob- 



scure and proves challenging to explain theoretically. The 
dispersion could be ascribed to a variability in the Li I 
line as a consequence of stellar activity, different mixing 
processes, presence or absence of circumstellar disks, bina- 
rity, or different rotation rates from star to star. Recently, 
Fernandez & Miranda ( 1998| ) have found that the Lil 
A6708A line in the WTT star V410Tau varies accord- 
ing to its rotational period. From Figs. [□] and [l(] we ob- 
serve that the region of the largest lithium scatter coin- 
cides with that of the strongest Ha emissions. This might 
indicate that som e hot continuum i s "ve iling" the optical 
spectra (Joy 1945|; Basri & Batalha 1990|; Basri, Martin & 



Bcrtout |1991|) , thereby affecting our pEW measurements. 
We note, however, that if any "veiling" exists around Ha 
and Li I in our spectra, it has to be small compared with 
that of many other CTT stars, because there is no clear 
correlation between strong Ha emission and low values of 
Lil pEWs (except for S Ori J053951.6-022248). There are 
other possible explanations for the significant Lil pEW 
scatter, such as different gravities (objects with low Lil 
pEWs might have lower gravities, and therefore, younger 
ages), and contamination by lithium-depleted interlopers. 



5. Discussion 

5.1. Radial velocity: binarity 

Lithium detections guarantee youth and the very likely 
membership of our sample in the a Orionis cluster. With 
radial velocities we may be able to study possible multi- 
plicity. However, the large uncertainties and having only 
one epoch of observations for the majority of the targets 
prevent us from carrying out a detailed analysis. In gen- 
eral, the radial velocities in Table || are in the interval 30- 
50kms _1 . Walter et al. ( |1998| ) obtained radial velocities of 
104 pre-main sequence stars within 30' from the a Orionis 
star. These authors find a sharp distribution peaking at 
around 25kms _1 and covering a range from lOkrns -1 up 
to 50 km s _1 . Our measurements are in full agreement with 
this wide radial velocity survey. Figure ^ depicts our ra- 
dial velocities against I magnitudes. Neither drift nor an 
increasing dispersion are obvious at the faintest magni- 
tudes, indicating that the very low mass stars and brown 
dwarfs of a Orionis are not yet affected by internal dynam- 
ical evolution, and that these objects still share the bulk 
motion of the group. 

Only the brown dwarf S Ori 45 clearly shows a rather 
discrepant radial velocity, which differs by more than 2.5 a 
with respect to the cluster mean velocity. With a mass es- 



timated at around 0.02 Mq (Bejar et al. |1999[ ), S Ori 45 
is the smallest object in our sample. It might belong to 
another kinematical group of young stars, like the Taurus 
star- forming region or the Gould Belt. On the basis of its 
multi- wavelength photometry and spectroscopy, S Ori 45 
is probably not a member of Taurus. The distance mod- 



CD 

CD 



o 




I magnitude 

Fig. 17. Radial velocities against / magnitudes. S Ori 45 
is not included in the figure (see text). 



ulus to Taurus is 5.76 (Wichmann et al. 1998), which 
would make S Ori 45 incredibly overluminous by 2.2 mag 
in the HR diagram. Guillout et al. ( [1998 ) and Alcala et 
al. ( 2000| ) have shown that the distribution of candidate 
members of the Gould Belt for the particular direction 
towards Orion lies at 200-300 pc from the Sun and well 
to the southwest of the Orion A cloud. This is relatively 
far away from a Orionis (>55pc). S Ori 45 fits the photo- 
metric and spectroscopic sequences of the a Orionis clus- 
ter very nicely (Bejar et al. 1999| , |2001 ), supporting its 
location in the Orion complex. Furthermore, this brown 
dwarf displays strong Ha emission and lithium in its at- 
mosphere, which is typical of ages much younger than that 
of the Gould Bel t (30- 80 Myr, Alcala et al. [2000| ; Moreno, 
Alfaro & Franco 1999| ), and it does not show a radial ve- 
locity consistent with membership in either Taurus or the 
Gould Belt. Alternatively, S Ori 45 might be a runaway 
object of the a Orionis cluster resulting from encounters 
with other cluster members; it may have been dynami- 
cally ejected from the multiple system where it originated 
(Kroupa |1998|; Portegies Zwart et al. |1999|; Reipurth & 



Clarke |200l| ; Boss [200 lQ , or S Ori 45 might be a brown 
dwarf close binary. So far none of these hypotheses can 
be discarded. Further radial velocity measurements are 
needed to assess the possible binary nature. If S Ori 45 is 
proved to be a spectroscopic binary, the dynamical masses 
of the components will be valuable for testing theoreti- 
cal evolutionary tracks at very young ages and substellar 
masses. 

From Fig. |k] we observe that r053820-0237 (M5) ap- 
pears remarkably overluminous with respect to the cluster 
photometric sequence. In addition, its radial velocity is the 
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Fig. 18. Surface curves for lithium depletions by factors 
of 3 (log AT (Li) = 2.5, upper panel) and 10 (log iV(Li) = 2.0, 
lower panel) as a function of age and mass. Models are 
taken from D'Antona & Mazzitclli ( |1994 , dashed line), 
D'Antona & Mazzitclli (1997, dotted line), Pinsonncault 
et al. ( |1990[ solid line) and Baraffe et al. ( |1998| , dash- 
dotted line). 



largest amongst our measurements. These two properties 
suggest that this star is an equal mass binary. 

5.2. The age of the a Orionis cluster 

We do not observe from Fig. |l6| that our a Orionis targets 
have undergone appreciable lithium destruction. Actually, 
the Li I curves of growth that neatly reproduce the ob- 
servations are those computed with the "initial" lithium 
abundance. More massive F- and G-type stars in the Orion 
complex have similar lithium contents (Cunha, Smith & 
Lambert 1995 ). The lower envelope to the distribution of 
the Li I pEWs shown in Fig. |l^ could be described by the 
log A^(Li) = 1.9 curve of growth, i.e., lithium depleted by 
about one order of magnitude. We shall discuss the likely 
age of the a Orionis cluster on the basis of no lithium 
destruction, and depletions by factors of 3 (logarithmic 
abundance of ~2.5dex) and 10 (logarithmic abundance of 
2.0 dex). 

According to various evolutionary models available in 
the literature, very low mass stars (M <0.3 Mq) burn 
lithium very efficiently by one order of magni tude a t ages 
older than 15Myr (D'A ntona & Maz zitclli |l994| , |l997 ; 



199J). Stars with 



Pinsonneault et al. |1990| Baraffe et al. 
masses in the interval 0.5-0.8 M Q do it in a shorter time 
scale. This is summarized in Fig. [l^, which shows sur- 
face curves for a given lithium abundance as a function of 



age and stellar mass. The age of the a Orionis cluster will 
be constrained by late-K and early-M stars. More massive 
members (M > 0.9 Mq) need longer times to deplete some 
lithium, so they are not useful for our purposes. 

Lithium depletion by a factor of 10 will impose a rather 
conservative upper limit on the age of the cluster. From 
Fig. ^ we infer that this upper limit is around lOMyr 



(based on Baraffe et al. 199S and Pinsonneault et al 



1990 models), because this is the time required by 0.6- 
0.8 Mq -stars to consume their lithium from initial abun- 
dance dow n to l og AT( Li) = 2.0. Models by D'Antona & 
Mazzitelli ( 1994 , 1997| ) predict values that are twice as 
young, i.e., around 5Myr. However, since no lithium de- 
pletion is apparent in any cluster member, it seems rea- 
sonable to establish shorter upper limits. If we adopt the 
surface curve corresponding to a factor of 3 lithium de- 
piction, the plausible oldest age o f the a Orionis cluster is 
8Myr (as g iven by Baraffe et al. 1998 and Pinsonneault 
et al. 1990] models). We have also inspected the lithium 
depletion tracks provided by ProfHtt & Michaud (1989) 
and Soderblom et al. ( |1998 ) obtaining very similar values. 
Our result fully agrees with the maximum age expected 
for the central, most massive cluster star to blow up as a 



supernova (Meynet et al. 1994). a Orionis low mass stars 
span an age range similar to that of the early-type mem- 
bers, i.e., the low and high mass populations are essen- 
tially coeval. Similar upper limits are found for other as- 
sociations in Orion, like the star-forming region around 
the A Orionis star (7-8 Myr, Mathieu, Dolan & Robert 
2001D , and Orio n WTT stars (Alcala, Chavarrfa-K., & 
Tcrranegra 1998D . We could adopt as the mean cluster 
age the oldest isochrone for which lithium is still preserved 
within 0.2 dex across the entire mass range. This occurs 
at roughly 2-4 Myr considering all models, a result in full 
consistency with previous analysis of theoretical isoch rone 
fitting to the observed photometry (Bejar et al. 199E). 



An additional constraint to the age of the cluster comes 
from the ratio of CTT stars to WTT stars. Based on 
strong Ha emission and the presence of forbidden emis- 
sion lines, this ratio turns out to be in the range 30- 
40% in a Orionis. Follow-up observations of our targets 
(mid-infrared, radio) are, however, desirable to confirm 
the presence of circumstellar disks. The ratio obtained 
in a Orionis is slightly smaller than that of younger re- 
gions, like the area around the Orion Molecular Cloud 



(ratio >40%, 1-3 Myr, Rebull et al. |2000J) , and consider- 
ably larger than the one of older clusters and associations, 
like the Sco-Cen OB association (ratio of 11%), whose pop- 
ulation of CTT stars, WTT stars and post-TTauri stars 
has been investigated by Martin Q1998 ). This author de- 
fines post -T Tauri stars as young, late- type stars that are 
burning lithium and display moderate Ha emission. The 
average age of the whole Sco-Cen OB association is in the 
range 5—15 Myr, as determined by de Geus, de Zeeuw &z 
Lub ( 1989 ). Wc do not find evidence for the existence of 
post-T Tauri stars in a Orionis, and hence, this cluster is 
essentially younger than the Sco-Cen OB association. 
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6. Summary and conclusions 



References 



We have presented intermediate- and low-resolution opti- 
cal spectra between 6100 A and 7000 A, covering Ha and 
Li I at A6708A, for a total of 25 low mass stars and 2 
brown dwarfs members of the a Orionis young star clus- 
ter. Spectral types have been derived and are found to 
be in the interval K6-M8.5, which corresponds to masses 
from 1.2 Mq down to 0.02 Mq after comparison with 



state-of-the-art evolu tionary models (Baraffe et al. |1998 ; 
Chabrier et al. |200C| ). We have measured radial velocities 
and pseudo-equivalent widths (pEWs) of the Ha and Li I 
atomic lines and find that all our targets show remarkable 
Ha emission and Lil A6708A in absorption. All radial ve- 
locities (except for one object) are consistent with mem- 
bership in the a Orionis cluster as well as in the Orion 
complex. The distribution of Ha and Lil pEWs against 
spectral type exhibits a large scatter at classes cooler than 
M3.5. This phenomenon occurs at the approximate mass, 
~0. 25 Mq, where low mass stars are expected to become 
fully convective. Some of our objects also show emissions 
of He I and forbidden emission lines of [O i] , [N n] and [S n] , 
probably indicating accretion from circumstellar disks. We 
infer that the likely rate of a Orionis low mass members 
resembling classical TTauri stars is in the range 30-40%, 
suggesting that the cluster is only a few Myr old. 

We note the intriguing case of the coolest object in our 
sample, SOri45, an M8.5-type brow n dw arf with a mass 
estimated at 0.02 Mq (Bejar et al. 1999 ). It has a very 
intense, variable Ha emission and lithium in absorption. 
Our tentative detection of forbidden emission lines of [N n] 
and [S n] suggests that S Ori 45 may have a cool, surround- 
ing disk from which it is accreting. This brown dwarf also 
displays a radial velocity that deviates significantly from 
the cluster mean velocity. 

We have also presented very recent computations of Li I 
A6708A curves of growth for low gravities (log g =4.0 and 
4.5), cool temperatures (T off = 4000-2600 K), and lithium 
abundances in the interval log ./V(Li) = 1.0-3.4. The dis- 
tribution of our observed Lil pEWs appears to be well 
reproduced by the theoretical pEWs computed for the cos- 
mic lithium abundance of logiVo(Li) — 3.1. This leads us 
to conclude that lithium has not yet been depleted in the 
a Orionis cluster. Therefore, after comparison to various 
lithium depletion curves available in the literature, we im- 
pose an upper limit to the cluster age of 8 Myr, while the 
most likely age is in the interval 2-4 Myr. 
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